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Abstract 

We consider quantum effects in a world with two coexisting symmetry phases, 
unbroken and spontaneously broken, as a result of a first order phase transition. The 
discrete symmetries of the problem are discussed in general. We compute the exact two- 
point Green function for a free fermion, when a thin wall separates the two phases. 
The Dirac propagator displays both massive and massless poles, and new CP-even 
phases resulting from the fermion reflection on the wall. We discuss the possible quark- 
antiquark CP asymmetries produced in the Standard Model(SM) for the academic 
T = case. General arguments indicate that an effect first appears at order aw in 
the reflection amplitude, as the wall acts as a source of momentum and the on-shell 
one-loop self-energy cannot be renormalized away. The asymmetries stem from the 
interference of the SM CP-odd couplings and the CP-even phases in the propagator. 
We perform a toy computation that indicates the type of GIM cancellations of the 
problem. The behaviour can be expressed in terms of two unitarity triangles. 
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1 Introduction 



In ref. [0 we presented recently a summary of our ideas and results on Standard Model (SM) 
baryogenesis, in the presence of a first order phase transition. The aim of the present work 
is to describe in detail our zero temperature analysis. The finite temperature T scenario is 



treated in the acompanying paper, ref. ||14|| . 



The baryon number to entropy ratio in the observed part of the universe is estimated to 
be tib/s ~ (4 — 6)10 _11 [[J. In 1967 A.D. Sakharov[[|] established the three building blocks 
required from any candidate theory of baryogenesis: 

(a) Baryon number violation, 

(b) C and CP violation, 

(c) Departure from thermal equilibrium. 

The Standard Model (SM) contains (a)0 and (b)||, while (c) could also be large enough 
1| , if a first order SU (2) x f/(l) phase transition took place in the evolution of the universe 
||. We will not enter the discussion on the latter: it will be assumed that a first order phase 
transition did take place, as strong as wished by the reader. An optimal sphaleron rate can 
be assumed as well. Our aim is to argue, on a quantitative estimation of the electroweak 
C and CP effects exclusively, that the current SM scenario is unable to explain the above 
mentioned baryon number to entropy ratio. 

Intuitive CP arguments lead to an asymmetry many orders of magnitude below observa- 
tion [|Kj[|rT|. Assume a total flux of baryonic current, where all quark flavours are equally 



weighted. An hypothetical CP- violating contribution in the SM with three generations^] 
should be proportional to 

si 2 s 2 s 3 c 1 c 2 c 3 ss {m 2 - m 2 )(m 2 - m 2 ){m 2 - m u 2 ) 

x(m b 2 - m 2 ){m b 2 - m 2 ){m 2 - m d 2 ), ( 1.1) 

times some power of the electroweak coupling constant, to be determined. In eq. ( | 1.1| ) the 
mixing angles and phase are the original Kobayashi-Maskawa ones||. In order to obtain a 
dimensionless result, the expression in ( | 1.1| ) has to be divided by the natural mass scale 
of the problem Mw or, at finite temperature, T, at the 12 th power. It follows|ll]] that the 
resulting asymmetry is negligible, %/s ~ 10~ 20 or smaller if other reduction factors are 
considered. When the physical boundary conditions are such that some specific flavour is 
singled out by nature (or experiment as, for instance, in the K — K system, electric dipole 
moment of the neutron, etc.), some of the fermionic mass differences in Q 1.1]) are no more 
compelling, and a stronger effect is possible. Nevertheless it seems difficult to compensate the 
10 orders of magnitude difference in the above reasoning. Such could be the case if the desired 
observable came about as a ratio of ( | 1.1|) with respect to another process which contained 
by itself some fermionic suppression factors: an example is the CP-violating e parameter in 
K decays, where the CP-violating amplitude is compared to a CP-conserving denominator, 
~Re (K — Kq), which is a GIM suppressed object. This is not what happens with the 
baryon asymmetry of the universe where, as shown below, although nature may single out 
some specific flavours, there is nothing to divide with (other than the overall normalization 
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to the total incoming flux, which bears no suppression factor). In this respect, the process 
has some analogy with the electric dipole moment of the neutron where the experiment 
preselects the up and down flavours, and of course the results are normalized to the total 
incoming flux of neutrons. 

On the light of the above considerations, the reader may be astonished that we give 
the problem further thought. Our motivation is that the study of quantum effects in the 
presence of a first order phase transition is rather new and delicate, and traditional intuition 
may fail. A SM explanation would be a very economical solution to the baryon asymmetry 
puzzle. To discard this possibility just on CP basis suggests that new physics is responsible 
for it, without any need to settle whether a first order phase transition is possible at all. 
Furthermore, the detailed solution in the SM gives the "know how" for addressing the issue 
in any theory beyond the standard one where CP violation is first present at the one-loop 
level. On top of the above, the authors of ref.[|l2|] have recently studied the issue in more 
detail and claim that, at finite temperature, the SM is close to produce enough CP violation 
as to explain the observed %/s ratio. 




Figure 1: Quarks scattering off a true vacuum bubble. Some notations in the text are sum- 
marised in the "zoomed-in" view. 

A first order phase transition can be described in terms of bubbles of "true" vacuum 
(with an inner vacuum expectation value of the Higgs field v 0) appearing and expanding 
in the preexisting "false" vacuum (with v = throughout). We can "zoom" into the vicinity 
of one of the bubbles, see Fig. [L| There, the curvature of its wall can be neglected and 
the world is divided into two zones: on the left hand side, say, v = 0; on the right v ^ 
and masses appear. The actual bubble expands from the broken phase (v ^ 0) towards 
the unbroken one (v = 0). We work in the wall rest frame in which the plasma flows in 
the opposite direction. Consider thus a baryonic flux hitting the wall from the unbroken 
phase. Far enough to the left no significant CP-violating effect is possible as all fermions are 
massless. In consequence, the heart of the problem lies in the reflection and transmission 
properties of quarks bumping on the bubble wall. CP violation distinguishes particles from 
antiparticles and it is a priori possible to obtain a CP asymmetry on the reflected baryonic 
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current, Aqp- The induced baryon asymmetry is at most np/s ~ 10 2 Aqp, in a very 



optimistic estimation of the non-CP ingredients [T3][l^ . 

The symmetries of the problem are analyzed in detail for a generic bubble. The analyt- 
ical results correspond to the thin wall scenario. The latter provides an adequate physical 
description for typical momentum of the incoming particles \p\ smaller than the inverse wall 



thickness I, i.e., \p\ <C l/l- For higher momenta the cutoff effects would show up[17], but it 



is reasonable to believe that the accuracy of the thin wall approximation produces an upper 
bound for the CP asymmetry. 

The precise questions to answer in the above framework are: 1) the nature of the physical 
process in terms of particles or quasi-particles responsible for CP violation, 2) the order in 
the electroweak coupling constant, aw, at which an effect first appears, 3) the dependence 
on the quark masses and the nature of the GIM cancellations involved. 

We shall consider the problem in two steps: T = in the present paper, and finite 



temperature case in the subsequent one [TJj]. The cosmological first order phase transition 
is a temperature effect. In order to disentangle the physical implications of the presence 
of a wall, with the consequent breaking of translation invariance, from the pure thermal 
effects, we consider here an hypothetical world at zero temperature but with two phases of 
spontaneous symmetry breaking. We will see that this academic model is illuminating. 

Intuition indicates that an existing CP violating effect already present at zero tempera- 
ture will diminish when the system is heated, because the effective v.e.v. of the Higgs field 
decreases and in consequence the fermion masses do as well (only the Yukawa couplings 
already present at T = remain unchanged). The expected decrease of CP asymmetries 
for increasing T follows, then, from the well known fact that the Kobayashi Maskawa CP 
violating effects of the standard model disappear when at least two fermion masses of the 
same charge vanish or become degenerate. This intuition can be misleading only if a new 
physical effect, absent at T = and relevant for the problem, appears at finite temperature. 
Treating first the T = case allows a clean analysis of the novel aspects of physics in a world 
with two phases of spontaneous symmetry breaking. The quantum mechanical problem is 
well defined and in particular the definition of particles, fields, in and out states, etc. is 
transparent. 

At T = the building blocks of the CP violating effect are threefold. First, the necessary 
CP-odd couplings of the Cabibbo-Kobayashi-Maskawa (CKM) matrix are at work. Second, 
there exist CP-even phases, equal for particles and antiparticles, which interfere with the 
pure CP-odd ones to make them observable. These are the reflection coefficients of a given 
particle hitting the wall from the unbroken phase. They are complex when the particle energy 
is smaller than its mass. Third, as argued below, the one loop self-energy of a particle in 
the presence of the wall cannot be completely renormalized away and results in physical 
transitions. The effect, thus, is present at order aw in amplitude. Such an effect is absent 
in the usual world with just one phase of spontaneous symmetry breaking, where the self- 
energy is renormalized away for an on-shell particle. The difference is easy to understand: 
the wall acts as an external source of momentum in the one-loop process. 

The truly and essential non-perturbative effect is the wall itself. The propagation of any 
particle of the SM spectrum should be exactly solved in its presence. And this we do for a 
free fermion, leading to a new Feynman propagator which replaces and generalizes the usual 
one. With this exact, non-perturbative tool, perturbation theory is then appropriate in the 
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gauge and Yukawa couplings of fermions to bosons, and the one loop computations can be 
performed. 

Strictly speaking, the gauge boson and Higgs particles propagators in the presence of 
the wall are required, and it is possible to compute them with a similar procedure |L6]. Fur- 



thermore, at T = the field theoretical problem is ill-defined in the unbroken phase where 
the fermions are massless, due to infrared divergences. Some cutoff is needed, although 
no dependence on it should remain when the cancellations between different diagrams con- 
tributing to a CP-violation effect are considered. By the time being, we work in a simplified 
case in which the one-loop effect itself is computed in one of the two phases. The final 
result describes a transition between any two flavours of the same charge, resulting in a CP 
violating baryonic flow for a given initial chirality, of order a\y 2 in the total rate. The GIM 
cancellations appear as powers of mf 6 /Myy 6 times logarithmic corrections, for the internal 
fermionic masses in loops. The dependence in the external masses is not trivial either, and 
its GIM implications are discussed in detail. It follows numerically that the CP factor is 
orders of magnitude smaller than the CP factor needed for baryon number generation. This 
T = model is academic, and we develope it to sharpen our tools for the physical finite 
temperature case [p~4]1 . 

The scope of the present paper goes beyond the particular issue of baryon number genera- 
tion in the SM. For instance, the above mentioned exact fermion propagator in the presence 
of a wall, which we derive, should be useful in other scenarios when a first order phase 
transition is present. The understanding at the quantum level of physical processes in the 
presence of phase transitions still is in its infancy, and we give it a modest try. CP violation 
in the SM is just an example of an effect which would disappear in a classical or semiclassical 
statistical physics approach. 



2 Notations and symmetries 

In this section we settle the formalism and express some general results. Consider several 
quark flavours in a gauge theory with a vacuum presenting a "wall" structure, parallel to the 
x— y plane. More precisely, assume that the Higgs field (or fields if there are several Higgses) 
has a vanishing v.e.v. for z — ► — oo and some constant v.e.v. for z — > oo. In between there 
is a "wall" with arbitrary width and profile. The interactions of quarks with the gauge 
and scalar bosons are given by some gauge invariant Lagrangian. As previously stated, the 
quark-boson interactions will be treated in perturbation, while the non-perturbative effect 
of the wall on the fermion propagation is solved exactly. 

As a first step let us consider the non-interacting case and treat in first quantization 
formalism the free Dirac Hamiltonian, in the presence of the wall. The second quantized 
fields in the "interaction representation" are built subsequently. 

The quarks/antiquarks hitting the wall from the unbroken phase are reflected or trans- 
mitted. The question is whether there exists a CP asymmetry that produces a different 
reflection probability for quarks and antiquarksi. 

3 Note that current conservation relates transmission to reflection, so that an asymmetry in the reflection 
probability automatically implies an asymmetry in the transmission probability. 
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In the unperturbed case, we will diagonalize the free Hamiltonian. The eigenstates pro- 
vide then the reflection and transmission coefficients. The fields and the quark propagators 
are made out of these eigenstates. Furthermore, they will be the building blocks of the in- 
coming (from the unbroken phase) and outgoing (back to the unbroken phase) wave packets 
when higher orders are considered later on. 

The symmetries of the problem and the conserved or partially conserved quantum num- 
bers are relevant to classify the eigenstates. It should be obvious that discrete symmetries 
are of particular relevance. In the second subsection we shall derive in the exact theory the 
consequences of the only unbroken discrete symmetry, CPT, and study in that context the 
effect of CP violation. Note that the consideration of a flux of particles approaching the wall 
with a non-zero average velocity (which will correspond later on to the physical process) 
breaks the CPT symmetry in a global sense, but the microscopic relations for two-particle 
transition amplitudes derived below will still hold. 

2.1 Free quarks in the presence of a wall. 

Let us consider a static "wall" parallel to the x — y plane. The Higgs vacuum expectation 
value depends only on the z coordinate, resulting in hermitian mass matrices 777(2:) and 7715(2:) 
for the quarks. The time independent Dirac Hamiltonian is then 

H = a ■ p + (3m(z) + if3^ 5 m 5 (z) , (2.1) 

where a« = 70 7« and (3 — 70. 

For the SM, in the basis where the mass matrices are diagonal, the study simplifies to the 
one flavour case and 777(2:), 7775(2:) become real. Furthermore, it is always possible to rotate 
777.5(2:) away through a chiral rotation. In other models, a diagonalization of both matrices 
is not always possible and the matrix 7775(2) may remain, for instance, in two-Higgs models 
when the relative phase between the corresponding v.e.v.'s changes with z (see for example 
ref[0). 

Let E,p x ,p y ,p z denote the four-momentum of the quark. For a static wall the particle 
energy, E, is conserved, while p z is not when the fermion crosses the wall or bounces back, due 
to lack of translational invariance. The system is symmetric under reflections with respect 
to the x — y plane and with respect to rotations around the z axis, the latter implying 
conservation of total angular momentum in the z direction, J z . It is also invariant under 
Lorentz boosts parallel to the x — y plane (K x and K y commute with H). From the two 
preceding symmetries follows a conserved quantum number for any wave function which is 
an eigenvector of E,p x and p y . Boost the reference frame to obtain p x = p y = 0. In that 
frame, the helicity states of the incoming plane waves in the unbroken phase correspond 
to eigenstates of J z = S z . This quantum number is more convenient to use than helicity 
because it is conserved. Denote generically by j z (j z = ±1/2) the corresponding eigenvalue 
of J z : it will be used when labeling the eigenstates of the Hamiltonian. 

Charge conjugation, C, is conserved when there exists a basis where both 77i(2:)and 7775(2:) 
are real matrices. 

Define P' as the product of parity times a rotation of angle tt around the y axis, i.e., 
parity with respect to the x — z plane, which will maintain the unbroken phase on the initial 
side of the z axis. P' is a symmetry when 7775 (z) = in some basis. CP' is conserved when a 
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basis exists such that m(z) + i r y 5 m^{z) is a real matrix. The unperturbed SM Hamiltonian 
is invariant under C, P' and CP'. In general, the following quantum numbers are conserved: 
p x ,p y , J Z ,E and CP'T. 

To be specific we have assumed that the unbroken phase is on the negative side of the z 
axis. We define the eigenstates ip n by 

E n iP n = Hij n ( 2.2) 

where the "energy" E n may be positive or negative, and n is a short hand notation for the 
quantum numbers that label the states. For the latter we choose n = p x ,p y ,j z ,E n ,f,a. 
p x ,Py and E n are conserved real quantum numbers. They can take any positive or negative 
value provided E\ > pi + p 2 y . j z is a discrete conserved quantum number defined above. 
/ is the flavour of the state in the broken region, a will be relevant when the dimension 
of the eigenspace is bigger than one, for fixed values of the other quantum numbers. The 
eigenspace has dimension one in the case of total reflection, i.e. when the energy is smaller 
than the height of the wall, \E n \ < m(+oo), where it is assumed that the basis is such that 
^5(00) = 0. In the opposite case, when the energy is larger than the height of the wall, it 
has dimension two. a serves then to label the state in any convenient basis. Examples will 
be given in the case of the thin wall. With the above tagging, two states with different labels 
are linearly independent. 

Given a complete set of orthonormal eigenstates of the Hamiltonian, fields are defined in 
the usual way, 

*(£, t) = E b n^n + (x)e-^+ 4 + J2 (x)e-^~\ ( 2.3) 

n+ n~ 

where n + (n~) label the positive (negative) energy states. b n + and d n + and their hermi- 
tian conjugates are annihilation and creation operators 4 with the usual anticommutation 
relations. The fields *f>(x,t) verify the Dirac equation in the presence of the wall. The 
completeness relation for the eigenstates ip n {x) forces upon them the standard Fermi-Dirac 
anticommutation relations. They are the Heisenberg quark fields in the presence of the wall, 
when the coupling of quarks to the gauge and scalar bosons are switched off: they constitute 
the starting point of usual perturbation theory. 

Strictly speaking, although we study transitions between two given fermions, the prop- 
agation of bosons in the presence of the wall is needed. The reason is that a CP-violating 
effect in the SM is at least a one-loop process. The bosonic propagators in the presence of 
the wall are thus pertinent and viable in analogy to the fermion treatment we are developing 
here[|I(J. For technical simplicity, we leave such a task for a forthcoming publication, as the 
subsequent modifications should not change neither the order in the electroweak coupling 
constant at which the effect first shows up, nor the type of GIM cancellations. 

Table [I] shows the transformation laws for the discrete symmetries of the system. 



2.2 Symmetries of the exact theories. 

Consider now the exact theory with the quark-boson interactions incorporated. The states 
in eq. (| 2.2 ) are no more eigenstates of the total Hamiltonian. The fields defined in eq. ( 



The annihilation operators verify 6 n |0 >= d n \0 >— 0, the state |0 > being the "vacuum" that contains 



the wall, i.e. with a Higgs expectation value that depends on z. 
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symmetry 


wave function 


field 


momentum 


helicity 


Identity 




#(z,t) 


P 


h 


P' 




i<jy 7o^(5, t) 


P 


-h 


C 




2 72 ^* (£, t) 


P 


h 


CP 1 




75** (x, t) 


P 


-h 


T 


~(J y i)* n {x) 


-a y V(x, -t) 


-P 


h 


TCP' 


-Vyl<o^n{x) 


-^^(x, -t) 


-P 


-h 



Table 1: Transformation laws for the wave functions and the fields under the discrete sym- 
metries of the wall, with p = (p x , —p y ,p z ) ■ In the case of the symmetries C, CP' and CP'T 
the transformed wave function corresponds to a negative energy state. However the trans- 
formed momenta and helicities we quote are those of the anti-particles, i.e. opposite to those 
for negative energy states. 

2.3) are not Heisenberg fields of the full theory, but they can be used as the quark fields 
in the "interaction representation". The discrete symmetries apply to the fields exactly as 
stated in Table [I]. In general all these symmetries are broken except TCP'. 

Consider, for definiteness, the unbroken phase. From Table [I] the TCP '-transformed of a 
state with flavour /, momentum p and helicity h is given by 

\f,p,h> TCP '=\f,-p,-h>, (2.4) 

where / corresponds to the antiquark with the same flavour. TCP' is an antiunitary opera- 
tion, and invariance under it implies that 



TCP' <b ^-iH(f 1 -t^ a> TCP> = (<6 | e -^( tl - t2 )| G>) * 

= <a|e- <H(t i-^ ) |6>, (2.5) 

where a, b are two given quark states and t[ = —t\ and t' 2 = —t%- The relation 
e _ifr (*i -t 2) = ( e -*#(*i-*2))t been used. Equations ( | 2.4| ) and ( | 2.5| ) imply, for instance 
for a top to charm transition, that 

< c, (p x , Py ,-p z ),R\e- lH ^~^\t,p,L >= 

<i,(-p x , Py ,-p t ),R \e' iH ^-^\c,(-p x ,p y ,p z ),L>, (2.6) 

and the analogous formula for all flavours and helicities. 

Letting t\ — > oo and t 2 — > — oo, two important consequences follow from this equation: 

• i) All CP' asymmetries obtained by summing over all flavours, helicities and incoming 
momenta will necessarily vanish fL5|: 



£ \A 2 (f, h,p mc -> /', -h,p out )\ 2 - \A(f, h,p mc -> f , -h,p out )\ 2 = 0, ( 2.7) 

f,f',h,P inc 

where A(f, h,p inc — > /', —h,p out ) represents the amplitude for an incoming particle / 
with helicity h and moemntum p mc to be reflected into an outgoing particle /' with 
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helicity —h and momentum p out = (p™ c , p™ c , —p l z nc , E). This strong result does not kill, 
however, the scenario for baryogenesis during the electroweak phase transition, since 
the latter is based on the role of sphalerons, which is totally different for left-handed 
and right-handed quarks. In the following we will consider CP' asymmetries summed 
over flavours and momenta but only for a given initial helicity. 

• ii) A straightforward consequence of equation | 2.7| is 

\A(f,L,p mc - f,R,p° ut )\ 2 - \A(f,R,p inc - f',L,p° u X = 

f,f',p inc 

£ \A(f, L,p™ - f',R,p out )\ 2 - \A(f,R,p inc - f',L,p° ut )\ 2 , ( 2.8) 

which is a useful simplification, as it allows to express a CP' asymmetry without any 
use of ant i- quarks. 

As we are interested in the breaking of CP' symmetry, it is useful to express the conse- 
quences of an hypothetical theory invariant under this transformation: 

< c, (p x ,p y ,-p z ,R\e- iH ^-^\t,p,L>= 

<c,(p x ,-p y ,-p z ,L\ e- iH ^^\i,(p x ,-p y , Pz ),R>, (2.9) 

and the analogous formulae for all flavours and helicities. 



3 The thin wall. 

From now on we consider the SM in the presence of a thin wall, which allows to perform 
easier analytic calculations. The unperturbed Hamiltonian is, for each flavour, 

H = a-p + (3m6{z). (3.1) 



3.1 Eigenstates 

Consider a positive energy fermion with given p x ,p y , E$ 

On the left of the wall (z < 0), an incoming particle has a z-momentum p z = \J E 2 — p 2 x — p^ > 
0, and helicity h. Define 

P inc = (p x ,Py,Pz,E). (3.2) 
Upon hitting the wall, a reflected particle has ^-momentum — p z and 

P out = (Px,Pv,-Pz,E), (3.3) 

while a transmitted one has ^-momentum p' z given by 



Pz = vpI ~~ m2 ^ pI > m2 



p' z = iJm? — p 2 z if p 2 z < m 2 . ( 3.4) 



6 E>Jp*+p*. 
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With these definitions, e tp ' z ' z is a falling exponential in the case of total reflection (p 2 z < 
m 2 ). We thus define the transmitted 4-momentum as 

P tr = (p x ,p y ,p' z ,E). (3.5) 

The above physical situation is described by the following eigenstate of the hamiltonian, 
which will be denoted "incoming" eigenstate, 

^(x) = (u h {p inc ) e ipinc - s + Ru h {p inc ) e ipmt A 9{-z) 

+ (1 + i?)u h (p inc )e^ tr ^ 6{z), h=j z . (3.6) 

It is directly related to an incoming wave packet. Uh(p mc ) is a solution of the Dirac equation 
in the unbroken phase, 

fi nc u h (p inc )=0, (3.7) 
and the reflection matrix R is given by 

R = -^t- ( 3.8) 

Pz + P' z 

The Dirac structure of R reflects the opposite chirality of the incoming and reflected 
state. Note as well that R is complex in the case of total reflection. Its imaginary phase 
does not change from particles to antiparticles, and constitutes an example of the CP-even 
phases which will later on interfere with the CP-odd ones to produce observable effects. 

It follows from the definition of R that {1+R)u\] ic {p mc ) is a solution of the Dirac equation 
in the broken phase, 

(f -m)(l + R)u h (p mc ) = , (3.9) 

as well as 

f ut Ru h {p mc ) = . (3.10) 

In fact, eq. ( | 3.9|) was used to compute R. In the case of total reflection, p 2 — m 2 < 0, 
the wave function ( | 3.6| ) is the only positive energy eigenstate. 

When transmission occurs, one more eigenstate is needed to span the eigenspace. A 
second linearly independent solution, to be called "outgoing" state (because directly related 
to the outgoing wave packets), is given by 

= (u h (p out ) e l P ° ut - s + R) u h (p out ) e % P inc A 6{-z) 

+ (l + R^)u h (p out )e l P br - s 6(z), h=-j z , (3.11) 

where 

P br = (p x ,p y ,-(P' Z )\E). (3.12) 

This wave function tp out is however not orthogonal to ip mc . An orthogonal eigenstate is 
provided by the solution coming from the broken phase: 



+ (1 + J)M s (p br )e^ 



^ out . 



3.13) 
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where J is the reflection matrix for a particle bouncing on the wall from the broken phase, 
given by 

1 + J= — + (3.14) 

Pz 

and s is a spin index dependent on j z , such that (1 + J)u s {p br ) = Uh{p out ) is a massless 
spinor with helicity h = —j z , and u s (p br ) satisfies the Dirac equation in the broken phase 

(/ r -mK(p br ) =0 . (3.15) 

For later reference, we will name the solution in eq. (| 3.13|) the "broken incoming" wave 



function. 

The negative energy solutions, ip- n , describe the propagation of antiparticles. In terms 

of the quantum numbers for positive energy antifermions, the corresponding wave functions 
are 

^™_ c (f) = (v h (p inc ) e l P mc -* + Rv h (p inc ) e'P ""^ 9(-z) 

+ (l + R)v h (p mc )e l P tr - s 6(z), h = Jz , (3.16) 

xf;™l(x) = (v h (p out ) e#° ui - 3 + I$v h (j? out ) (?P in °A 9{-z) 

+(1 + i? f ) v h (p out ) br ' s 6(z), h = -j z , ( 3.17) 



and 



(&) = \IPz/p' z (v s (p br )e ij; br ' 3 + Jv s (p bT )e ii! tr '^0{z) 
+ (1 + J)v s (p br )e i P ° ut - s 6(z)], 



3.18) 



where R = —R and (1 + J) = p' z /p z (l + R). Vh(p mc ) is a negative energy Dirac spinor 
verifying the equation 

f nc v h (p inc ) = 0. (3.19) 

with the helicity h defined as a ■ pvh(p) = —hvh{p)- The spin s in eq. (| 3.18|) is defined 
as for eq. ( | 3 . 1 3| ) , in such a way that (1 + J)v s (p br ) = Vh{p out ) with helicity h = —j z . It 
follows from the definition of R that (1 + R)vh(p mc ) is a solution of the Dirac equation in 
the broken phase: 

U> tr + m)(l + R)v h (p inc )=Q, (3.20) 

as well as 

f ut Rv h (p mc ) = 0, (3-21) 

and 

(f r + m)v s (p br ) =0. (3.22) 
It is convenient to consider the Fourier transformed wave functions 

Hq) = J ^e^(z). ( 3.23) 
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From now on it will be assumed that the boost leading to p x = p y = has been performed. 
We will only write here the z dependence, since the x — y part is trivial (S(q x )5(q y )). The 
results are 

= I + fe) - ^TT, - -(J - ti - *) ) Uh{rU% (3 - 25) 

= Up-^T) - ,TT^ - ^5^}^' ( 3 - 27) 

= 2^ U - L + ie) - - -(J - - *) ^ (3 - 28) 

^-fe) = ^V / ^M(^^ ; / ^-t W^), (3.29) 

2ztt v I -p' z - (q z - ie) p' z - (q z - ie) -p z - (q z + ie) J 

with h = +j z in ?p™± and h = — j 2 in ip°±. 
3.2 Wave packets 

The physical process under study is described by a wave packet coming from the unbroken 
phase, bouncing on the wall, and generating reflected and transmitted wave packets. The 
transmission probability follows from the reflection one. It is thus sufficient to concentrate 
on reflection properties. 

Each wave packet has fixed flavour and helicity h (equivalent in the unbroken phase to 
jzi —jzi respectively for the incoming, outgoing wave packets). Its energy and momenta 
are clustered around average values. Somewhere in the process the electroweak interactions 
will act as a perturbation, and induce flavour changes through loops, as will be discussed in 
Section [5]. 

The x and y components of the wave packet are not modified by the wall, and will be 
ignored in what follows. Consider the following incoming wave packet, approaching the wall 
from the unbroken phase, 

P(p mc ,z,t) = N J dk z e-^- pTC)2 ^ e ik ^ z+z -^ r) u h {k) 

„ e i P ^(z+Z-t-T) e - (z+Z ^ T)€ Uh tf™^ (3.30) 

where N is a normalization constant, d denotes the spatial extension of the wave packet, 
which is located at time ~ — T (T > /) around position ~ — Z [Z > 0), assuming Z/d » 1, 
T/d » 1. d is chosen so that p l z nc d » 1. Terms exponentially suppressed by exp [— (p z nc d) 2 } 
have been neglected in the second line of eq. ( | 3.30| ). 

At t ~ — T, the wave packet is an almost monochromatic wave located far from the wall 
in the unbroken phase, with group velocity pointing towards the wall. To study its time 
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evolution when approaching the wall, an expansion on the eigenstates defined in subsection 
3~T| is convenient. Up to exponentially suppressed terms, the result is 



P(p inc , z,t)=N f dk z e -^-^" c ) 2 # i(t%f(z) e~ ik ^ t+r) . ( 3.31) 

In other words, the incoming wave packet totally expands on the eigenstates denoted by 
ip l n+i thus justifying their name. Analogous conclusions hold for the outgoing wave packet 
in the unbroken phase. It follows that the physical amplitude can be expressed at order aw 
as 

A(i -> /) =< 1%*\H M \rl>*r >, ( 3.32) 

where H int is the interaction Hamiltonian induced by electroweak loops. 



4 The quark propagator. 

The free quark propagator is defined as usual: 

S(x,y)=<0\T(*(x)*(y))\0>. (4.1) 

It is not translational invariant due to the presence of the wall, and verifies the equation 

(d xo + iH{x))S(x, y) = 7 (5 4 (x -y), ( 4.2) 

where H(x) is the Hamiltonian in eq. ( | 3.1| ). Furthermore, the time ordered product in 
equation (| 4.1| ) insures the Feynman boundary conditions. From eq. ( | 2.3|) , we get 

S(x,y) l0 = (6(x - y ) ^ iP n+ (x)-iPl + ($)e- iE ^-v°) - 

n+ 

8{Vo ~ x ) ^ i>nAZ)^n-(y)e- iEn - {x °- yo) ), ( 4-3) 

n~ 

where E n ± is the positive/negative energy of the eigenstate, cf. eq. ( | 2.2|) , and where the 
definition 

i = E (y) J d Vxdv y dE n ± ( 4.4) 

n± 3z,f,a 

has been used. 

The Fourier transformed propagator depends on two momenta, q l and q*, unlike in the 
translationally invariant case, 

q f = (E,q x ,q y ,qf), 

q i = (E,q x ,q y ,qi). (4.5) 



We can write 



x 



e 



iqi^+i^-iq^-i^-iqy^'y-^+iE^-^o) ^ / ^ 
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where translation invariance in x, y and time directions has been used. From now on we use 
the liberty of boosting in the x, y plane to bring q x and q y to 0. At the end of the section 
we will return to the general case. 

In terms of the momentum-space wave functions, the above translates into 



-1 



%Who = — T w 



c + c (^) (fete: 

+fe c te/) (fe c (?. 

+fe+te 2 7 ) (fete 

+$5r-(rf) (fete 



-1 


E- 


-E n+ +ie 




-1 


E- 


-E _-ie 




n 




-l 



E-E 

n 

-l 



E-E _-t 



4.7) 



Notice in the above expression that an orthonormal basis for the wave functions was used, 
i.e., the one spanned by ijj mc and ^) hr ', given in eqs. ( | 3.24Q , ( | 3.27|) , (| 3.26|) and (| 3.29|) . In 
particular, orthogonality requires to restrict the energy integration on the broken phase to 
E > m. We have explicitly verified the completeness relation for our system of eigenstatesi, 



fe+teD (fe c te:)) f + fe+tef) (fe + te:) x1 



fe-tef) (fe c teD) f + fetef) (fete:)) 1 



(4.8) 



(2tt): 



The obvious consequence of eqs. (| 2.2|) and (| 4.8|) is the Green's function equation for the 
propagator: 

- i(E - H) Stf, q l ho = (2n) 3 5(q{ - £). ( 4.9) 

The Fourier transform of eq. ( | 4.9| ) is: 



4.10) 



Details on the computation of the propagator can be found in Appendix [A]. Let us 
summarise the results in 1+1 dimensions (time and ^-directions). An interesting interme- 
diate step is given by the expressions which follow. The contribution of states propagating 
exclusively on the left-hand side of the wall (v = 0, "unbroken" region) is 



S(q f ,q l )iefao = {— f 

lai. — a 



q{ - q\ + ie (sE + ie) 2 - (qi + ie) 2 
1 1 s(l + sa 



(E + q f z a z 



sE + ie) - (q f z + ie) sE-qi + 



ie 



+ 



1 



'1 - sa. 



smjo 



sE + ql + ie sE-q\ + ie 



sE + p' 



4.11) 



6 The negative energy eigenfunctions describe the propagation of antiparticles: in the notation presently 
used an initial antiparticle will have momentum q* , while a final antiparticle has momentum q i . 
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The contribution from states propagating exclusively on the right-hand side of the wall (v.e.v. 
7^ 0), that is, in the "broken" region is 

S{q f , q'hghtlo = {— } , , 2 f (E + q{a z + m 7o ) 

L qi - q z ~ ie \Pz + lt ) - [—Qz + lt ) 

1 1 Elf p' m \ 

1 1 ^ z -a z + —7o 



(p> z + iE) + (ql -ie)p' z + ql + iep' z 2\ E z E 



p> z -ql + iep' z + ql + tep' z 2\ E z E J sE + p' 
Finally, the contribution from states propagating across the wall is 

S{q ; , q jacrosslO = S< - 7 ; — — 1 + , 

yp' z - qi + ie sE - q\ + ie \ sE + p' 0/ 

i i (!-*»,) r 1+ . . (413) 



sS + g^ + ie^ + 9j + * e 2 V sE + p' 
where s is a parameter which takes values +1 or —1 for positive or negative energy particles, 
respectively. p' Q is defined as follows 

p' Q = +y/E 2 -m 2 + ie, (4.14) 



and the factor srwyo/isE + p' ) originates from the reflection matrix R, eq. (| 3.8|) . 



Eqs. (| 4.11] )-( |4.13 ) show sometimes two ie terms with opposite signs present in the same 



denominator, for instance l/((sE + ie) — (q{ + ie)) in the second line of eq. (| 4. lip . The 



reason for this phenomenon is clear. The first ie in our example is there to implement the 
Feynman boundary conditions for the propagator in time direction. The second ie specifies 
in which spatial phase (z < or z > 0) this terms acts. This may seem ambiguous: is 
the pole located in the upper or lower half plane? Both prescription are not equivalent, 
unless the residue at the pole vanishes. In other words, since l/(x ± ie) = P(l/x) =F tt8(x), 
the ambiguity disappears only if the factor multiplying 5(x) vanishes at x — 0. A close 
scrutiny of eqs. (| 4.11|) - (| 4.13|) shows that this is indeed the case. For instance, the residue 
at sE = q{ in the second line of eq. ( | 4. 1 1|) is exactly compensated by the residue of the 



equally "ambiguous" term in the first line of eq. ( | 4.11Q . It results that the sum of both 
terms is not ambiguous. It is possible to chose in a consistent way any sign for the ie in the 
first and second lines of eq. ( | 4.11]) . Flipping from one convention to the other exchanges a 
term proportional to S(sE — ql) between the first and second line of eq. ( |4.11| ), and the sum 
does not change. For later use denote "time ie" the usual Feynman convention and "space 
ie" convention the other one. In the following the "time ie" convention will be assumed 
unless specified otherwise. 



The propagator is given by the sum of eqs. (| 4. lip , (| 4.12 ) and (| 4.13 ) above 



S'teVho = [S(q f ,q%ft + S(q f ,q% ight + S (q^ q') across ) 7o . ( 4.15) 



In Appendix |A] it is explicitly shown that eq. ( 4.15p verifies the Dirac equation. Using the 



relations ( | 5.38| - |~5.39[ ), as well as the above mentioned freedom to shift certain poles, it is 
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possible to combine the various terms in a more compact and familiar way. One can either 
show that 



1 

¥ 
i 



i 



i 



qt - q l z + ie f qi - q\ - ie f ~ m 

m 



+ 



"1 + sa z ~ 




ms7o 


2 




E + p' z _ 



— m 



ms7o 
E + p' z 



1 — sa. 



— m) 
m 



«7o- 



— m) 



4.16) 



which simplifies the action of the g*-Dirac operator, or that 

r 1 1 1 



S(q,q 



ql - ql + ie $ f q{ — q\ — ie\ 



m 



m 



i f (i f -m) 



ms7o 
E + p' z 



1 — sa, 



ni 



m) 



1 + sa z 



ms7o 
E + p' z 



s7o- 



m 



4.17) 



which prepares for the action of the g-^-Dirac operator. In fact, eq. ( | 4.16| ) is the CP' 
conjugate of ( | 4.17| ). An elegant expression, because more symmetric in the q l and q* 
dependences, can be obtained combining ( | 4.16| ) and ( | 4.17| ) and is given by 



S(q f ,ql = ~l/2{ - f 1 (■ 
^qi - ql + ie \\ 



1 1 



ql -q\- ie 



1 1 \ 1111 
+ I + — 7^77 - -Jllz^r 



m 



i f (i f -m) 



_ ms7o 
E + P'} 



m i % $ fa 1 — m 



sa z ) 



m 



s7o- 



m) 



( 4.18) 



In ( |4.16| ), ( |4.17| ) and ( |4.18| ): all ^'s in the denominators are understood as regularized "a 
la Feynman", i.e. as q" + ie, the "time ie" convention. This convention allows to check easily 
that our propagator obeys Feynman boundary conditions, i.e. that the positive frequencies 
(q > 0) are "retarded", and the negative ones (q < 0) "advanced". 

If the "space ie" convention was used instead, then the regularisation implies to consider 



all massless poles (1 



with q{. — > q{. + ie and q l z — > q % z —ie. Upon Fourier transform it 



corresponds to 9(—Zi) and 6{—Zf) respectively, i.e. it locates the massless poles in the unbro- 
ken region. Analogously, all massive poles (l/(^ — rrif), l/(jf — rrii)) are to be understood 
as ql — > ql — ie and q l z — > q l z + ie, locating them in the broken region. In this convention 
it is easy to check that the propagator in eqs. ( | 4.16Q , ( | 4.17| ) and 
equation (| 4.10Q . 



4.181) verifies the Dirac 
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For some terms in eqs. ( 4.16| ), ( | 4.17|) and ( 4.18|) the "time ie" and the "space ie" 
conventions are apparently in conflict. For instance, 



+ 2i^6{q{)5{ql - (g. 



f\ 2 \ 



4.19) 



The two conventions differ by the last term on this equation. However, this term, as all 
matching terms of this type, do cancel among the different factors in any of the eqs. ( |4.16 ), 
( [4171 ) and ( |TT8l) . This cancellat ion is in fact the same phenomenon as the cancellation of 
the residues corresponding to "ambiguous ie's" that we discussed after eqs. (| 4.11|) , ( | 4.12| ) 
and (| 4.13| ). We are thus lead to the striking conclusion that the "time" and the "space" ie 
conventions are strictly equivalent in eqs. ([JTT^A ( 1 4-lV an d ( 1 

This equivalence is physically understandable. It means that all the poles for which the 
two conventions give conflicting "ie" signs have a vanishing residue: they correspond to 
physically non acceptable asymptotic states. For instance, consider a massless pole with a 
positive q{ in eq. ( | A.7| ). Massless asymptotic states can only exist in the unbroken phase, 
i.e. on the left side of the wall. A final state in the unbroken phase must necessarily move 
away from the wall, i.e. with a negative q{. It is then physically welcome that the residue 
vanishes when q{ is positive. All other apparent conflicts are resolved in a similar way. 

We shall denote the first line of expressions (| 4.16| ), (| 4.17|) and (| 4.18| ) as the "non- 
homogeneous" part. The rest of the propagator in any of its versions will be denoted "ho- 
mogeneous part" , as the latter is a solution of the Dirac equation without a source, contrary 
to the former. By the same token, it is possible to verify that the advanced and retarded 
Green functions in the presence of a wall differ from those for the usual Dirac propagator by 
terms which vanish upon the action of the Dirac operator. 

The presence of both massless and massive poles reflects the fermion sailing between the 
unbroken and broken phases. It is trivial to verify in eqs. ( | 4.16|) , ( | 4.17|) and (| 4.18|) that 
the usual ijj) Feynman propagator for a massless fermion is recovered. Notice as well that, 
besides the usual ie dependence which corresponds to absorptive contributions for on-shell 
fermions, new imaginary phases appear which are generated by the reflection coefficients of 
the fermion when E 2 < m 2 . This CP-even phases, present for on-shell and off-shell fermions, 
are induced by the reflection matrix, eq.( | 3.8| ). 

For completeness we give as well the fermion propagator in the presence of the wall for 
the 4-dimensional case, i.e., when q x and q y have not been boosted to zero values: 



S{qf,J) = -l/2{ f 1 ( 
<qt — q\ + ie \ 




m 



m) 



- m g" qj q' — m 
E~fo ~ Qxlx ~ q y ly + ^gr(£7o - Qxlx ~ q y l y + s^ z ) 

Pz 



rn 



— m) 



4.20) 



where p z = , E 2 — q 2 — q 2 + ie and p' z = w E 2 — q 2 — q 2 — m 2 + ie. 
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During the (long) preparation of this manuscript, we became aware of ref.[nj which 
derives the Dirac propagator in presence of a thick wall and uses it to compute the CP- 
asymmetry in the reflection process. The source of CP-violation they use is not a SM loop 
but rather the tree-level axial quark mass in a 2 scalar model. Although this difference makes 
the comparison between both works delicate, their results seem to support our statements 
that increasing the wall thickness would further reduce the asymmetry. 



5 One loop computation. 

The propagators in the presence of the wall will be depicted by a line screened by crosses (see 
Fig. HD, while the usual ones are drawn as solid lines. Consider a one- loop transition between 
two external flavours, i and and /, such as the ones in Fig. |^. The analogous diagrams with 
full lines, i.e., with the usual Feynman propagators, are completely renormalized away for 
on-shell external states. This will be changed by the presence of the wall, which acts on the 
diagram as a external source of momentum. 



W 



H + ,H C 



(a) 



(b) 



(c) 




Figure 2: The kinky propagators for quarks, W's and scalars can be assembled into loop 
diagrams. Vertices, being local, are unchanged in this formalism. 

Let us review the properties expected in general, independently of the details of a given 
calculation. The process under consideration describes an asymptotically massless quark i, 
approaching the wall from the unbroken phase, which is reflected into the same phase as an 
asymptotic massless / quark. A flip of chirality has thus necessarily occurred. The flavour 
change happens through one electroweak loop, and we will denote by M the internal fermion 
mass. 

First of all, the amplitude should not vanish at order aw There is no symmetry reason 
known to us implying that the coefficient of the CP-violating operator should be zero at this 
order. 

Second, the amplitude should vanish at least as miUif when the internal loop is computed 
in just one of the two phases, and Lorentz covariance is preserved. A priori this behaviour is 
not mandatory, as just one chirality flip is needed, and an amplitude vanishing with masses 
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as either or would seem adequate. But there are more constraints: in the unbroken phase, 
the computation of the internal loop should behave as 



(5.1) 

due to Lorentz covariance. This loop is to be inserted on the propagator in the presence 
of the wall. Imagine that reflection has happened before the loop insertion (providing an 
mi factor), but not afterwards. The outgoing state is then a plane wave verifying the Dirac 
equation for a massless particle. The action of eq.( | 5.1| ) on it gives zero. In the opposite 
case, when reflection happens only after the loop insertion, the analogous argument for the 
incoming massless state holds as well. Consequently, the action of the wall must be present 
before and after the internal loop insertion, and a non trivial dependence of the amplitude 
on both rrii and mj results. The total effect should go to zero as some positive power of 
both of them, with an odd overall dependence on external masses, since chirality flips upon 
reflection. The argument can be generalized to the situation where the broken phase is 
considered inside the loop, the difference with eq. ( | 5.1| ) being terms independent of but 
proportional to the external masses, as it will be exemplified below. 

Third, it is possible to argue the type of GIM cancellations for the quark masses inside the 
loop. The relevant terms in the reflection amplitude A(i — > /) correspond to the interference 
of diagrams with two different internal quark masses, M and M', to be summed upon. Each 
individual diagram can be written as@ 



M 2 MM 4 M 
A(i - /) oc F(M) =a + — r (b + dog-—) + J r r (d + dog-—) + , ( 5.2) 

1V1 W Mw M w Mw 



with a, b, c, d,e... depending only on external masses and momenta. Due to CKM unitarity 
(see eq. ( | 5.25| )), the CP observable is proportional to J2m,m' Im[F(M)F(M')*]. As each 



term in this sum is an antisymmetric function of M and M', the squared terms in the devel- 
opment (| 5. 2D of the product cancel, while the cross-terms involving the constant piece a fall 



out once the sum is performed. The lowest order contribution is thus ~ J2 M lu 2 log(M 2 ) 

2 / 

In the practical example below, we find M . M log M' 2 , i.e. a further internal mass suppres- 

M w 

sion. It is important to notice that, whenever only the unbroken phase is considered inside 
the electroweak loop, c = d = e = 0, because the fermionic mass dependence stems from 
pure Yukawa couplings. No antisymmetric function in M, M' is viable, and the effect should 
vanish at order in total rate. There is no reason, though, to neglect the broken phase, 
and we expect an O(a^) contribution. 

Finally, the observable should be proportional to the interference of the imaginary parts 
of the CKM couplings with the CP-even phases produced in the reflection on the wall. 

The above considerations apply as well for the thick wall scenario. Some remarks specif- 
ically related to the wall thickness are appropriate, before entering into the details of our 
computation. The internal electroweak loop is in general a complicate function of q 2 , and 
thus non-local. Its typical "inverse-size" is of the order of the dimensional parameters in- 
volved which means > Mw, as we will see. To neglect its non-local character is only adequate 



7 This dependence is in general a complicate function. For transparency, only the behaviour for M < M\y 
is described here. As the internal loop is IR convergent when M — > 0, no pure logM terms can appear. 
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for wall thicknesses / larger than the loop "size", I » M^} 1 . For a thin wall, I — > 0, a lo- 
cal approximation for the internal loop is incorrect. Arbitrary large particle momenta are 
relevant, and it will be shown that indeed the non-local character of the internal loop is 
important. The role played by the latter effects in our calculation suggests that an even 
smaller result would follow in a more realistic thick wall computation, I >> , where a 
local internal loop could be a good approximation. 

We will see that the above four general characteristics do appear in the toy computation 
given below. 



5.1 One particle irreducible loop 

The complete computation of the diagrams in Fig. ^| requires not only the knowledge of 
the fermion propagator in the presence of the wall, but the corresponding one for gauge and 
Higgs bosons as well. They can be obtained in complete analogy with the computations of 
the previous section. However, they need rather lengthy calculations. Our goal here is rather 
to get a feeling of how the different building blocks work together than to solve exactly a 
problem which anyhow is academic to a large extent. We have chosen to simplify our task, 
leaving the full calculation for a forthcoming publication. What we compute in fact are the 
diagrams in Fig. [| with the internal loop computed in the broken phase. Notice that a 
similar choice has been made in recent publications which try to solve the problem at finite 



temperature [12], where only the unbroken phase was considered inside the electroweak loops. 




Figure 3: The subset of Feynmann diagrams used in this paper. 



Let us then consider the broken phase inside the loop. In the 't Hooft Feynman gauge 
the integral to perform can be written in general as 

= c t di \{ q -ky- a M^){k*-M?y ( 5 - 3) 

where g M is the external 4-momentum and a, b and c are independent. A indicates that the 
integral has been regularized someway, for instance by dimensional regularization. 

We choose to substract at q 2 = 0. Because Lorentz covariance is preserved at T = 0, the 
result of the substraction will have the general form 
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where A and B are matrices with 1 and/or 75 components, A subst (q 2 ) = A(q 2 ) — A(0) and 
B subst -(q 2 ) = B(q 2 ) - B(0). Finally, 

/-• (*.) = A jf dX /; Mia, + b) ^-0L^L__ . ( 5. 5) 

Notice that with, the above choice of substraction, the internal loop is completely renor- 
malized away for an on-shell massless fermion which suffers no interactions with the wall. 
Concentrate for the moment in the case i, f = d, s or b. The scalar-exchange contribution is 
then 

9 2 71-2 t^* ts f 1 j f 1 j ^M 1 - x ) 



2M^(2vr) 4 " Jo Jo y q 2 xy{l - x) - M^x - Mf(l - x) 
{x^niinifR + MfL) - M 2 (m f L + m^)}, 



5.6) 



where L and R denote the chiral projectors, L — (1 — 75)/% R = (1 + j 5 )/2. K describes the 
CKM couplings of the quarks and Mi is the mass of the internal quark of Z-th generation. 
The W-exchange contribution to the substracted internal loop is 

7r 2 r _ T ^ f 1 f 1 q 2 x(l - x) 

5.7) 



I$ bst -(q„) = g 2 -—K* u K u I dx / dy- — ^ ^ -x 

w y (2tt) 4 !/ Jo Jo q 2 xy(l-x)-M^x-Mf(l-x) 



The above represent only one contribution among many that should be considered in 
the full calculation. When this internal loop is inserted below on our propagator, gauge 
invariance is preserved, as the fermion interaction with the wall is of Yukawa type and thus 
SU(3) x SU(2) x U(l) gauge invariant. In the full theory, the substraction procedure will 
be more complicated, though!. Our aim in a toy computation is to get an idea of how the 
ingredients of CP vilolation, GIM mechanism, and wall reflection combine together, and to 
settle useful calculational tools for this type of problem. 

5.2 Loop insertion in the propagator 

The diagrams in Fig. |3] are computed inserting the internal loops given above on the fermion 
propagator in the presence of the wall, derived in the preceding section. The desired am- 
plitude should describe the transition of an incoming massless flavour i into an outgoing 
massless flavour /. The selection of these asymptotic states is obtained projecting onto the 
massless poles of the propagator, for the 4-momenta appearing at both extremes of Fig. S. 



8 What has been done above is close to the complete result in a limited extension of the standard model, 
with two Higgs's, one coupled to down quarks and the other (<3? u ) to up quarks. We assume the following 
expectation values: < $<j(z) >= v e i9(z) and < ^ u {z) >— v u , with Vd <C v u . In this model the equation ( 



5.7) turns out to become the exact contribution from W exchange, while eq. ( 5.6 ), or rather its Fourier 
transform, is modified by the following substitution: m., — * mi6(zi), rrif — > mfO(zf), where Zi (zf) is the 
coordinate of the incoming (outgoing) vertex of the loop. An analysis of the changes in the computations 
runs parallel to the one performed above, but it does not seem pertinent to pursue this model here. 
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The relevant contributions of Fig.(|3|) are then those terms which contain both poles in 
l/</f and l/jf^. It can be shown that the non-homogeneous part of the fermion propagator 
in the presence of the wall, gives no contribution of this type. That is, among the four 
possible combinations of homogeneous and non homogeneous on both sides of the internal 
loop, only the homogeneous-internal loop-homogeneous contributes, yielding for a positive 
energy particle: 



A{i - /) 



r dq z 1 


"1 + a z ~ 
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E + p'J 
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where 
and 



Pz 
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(5.9) 
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p z " = +yE 2 — m 2 f + ie. 

The q z integration reflects the arbitrary off-shellness of the quarks upon reflection on the 
wall. The relevant integrals are the following ones: 

dq z 1 1 1 



and 
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where 
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( 5-14) 

After some lenghty algebraic manipulations, the resulting amplitude for an incoming left 
asymptotic quark can be written as 
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5.17) 
5.18) 
5.19) 



Notice that Qo> obtained from the convolution of the electroweak loop with the free 
fermion propagator in the presence of the wall, is in general complex, and will be one of 
the relevant CP-even phases which will interfere with the CKM complex couplings. This 
new phase stems from the non-local character of the internal electroweak loop, and would 
have disappeared in a linear approximation. Other contributions of the same type should 
be present in a full calculation with the wall accounted for inside the loop, and although 
cancellations are possible, they are not mandatory by any symmetry argument. 

For an incident left-handed particle, the result is obtained from eqs. (| 5.15 ) to (| 5.19 ) 
replacing everywhere by mj and viceversa (and consequently p' z l by p'J). The CKM 
couplings remain identical. 

Consider transitions between charge —1/3 quarks. Concentrate as well on quark energies 
below or equal to Mw- The approximation m 2 ,mj,E 2 < Mf is then justified, and eq. ( 
5.15| ) simplifies to 
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where a, b, c, d, are given by eq. ( | 5.19Q , neglecting terms in m 2 / Mfy ,m 2 / Mfy . 
The amplitude can be rewritten as 



A{i R -> f L ) 
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2 (2tt) 4 

The function F(M) is defined as follows: 

1 



is ts* 1 milll S t?i tr\ 1 + a * T 
K ii K if-u^7T- F ( M ) = — L - 
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and the integrals Ij.,2,3,4 are 



I n (M) = 5 n / dx 



x(l — x) 



X + M^fi X ^ 



5.24) 



with 5 1 = 3, 5 3 = M 2 /M^, 5 2 = 35 4 /2 = 2 + M 2 /M^, (3 n = (0,1,0,1) and ln = 
(1,1,3/2,3/2). 

Consider i 7^ /, i.e., a flavour changing transition. The CP- violating observable is pro- 
portional to 



Y\m R - h)\ 2 - \A(j L - 1 
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where A(i — > /) is obtained from A(i — > /) taking the residues of 1/jf and 1/^ at g — 
= —5/. A(i — > /) corresponds to the matrix element of the effective Hamiltonian between 

states \i > and < /|, in a normalization that will be specified in the next subsection. 

Eq. ( I 5.25Q shows explicitly the interference between the CP-odd phases of the CKM 

couplings K and the CP-even phases of F. The latter can have two different origins: either 

p' z when i or / are totally reflected, or Q , which leaves a trace of the non- locality of the 

internal loop in the i factor before I 3 and I4. 
It is well known that 
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where 77 = 1 if/ — i = V — I (mod 3), 77 = —1 when / — % ^ V — I (mod 3) and where J 
is twice the area of the CKM unitarity triangle, J = CiC2C 3 s 2 s 2 s 3 ss. Let us now consider in 
the complex plane the triangle spanned by the three numbers F(Mi),l = u,c,t. Let B be 
the area of this triangle with a sign +1 (-1) for a cyclic clockwise (anti-clockwise) ordering 
of u, c, t around the triangle. Then 
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This expression exhibits clearly the GIM mechanism for the internal flavours, since the 
area vanishes obviously whenever two points coincide. Additional vanishing happens when 
three points are aligned. The function F depends both on the internal and external masses 
as can be seen from eq. ( | 5.23| ). In order to exhibit the GIM mechanism for the external 



masses as well, Im 



F{Mi)F*{M v ) 



can be expanded, leading to: 
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In this expression, the function S derives from the integrals /i, 2,3,4 m e Q- ( 5.24 ), and contains 
the dependence on the masses for the quarks inside the loop, 



S jk (M, M') = I 3 {M)I k {M>) - Ij(M')I k (M), ( 5.29) 

and the bjk are antisymmetric functions in j, k. They contain the dependence on the masses 
of the incoming and outgoing quarks. In the approximation used from eq. ( | 5.20| ) on, they 
are equal to 

b 12 (E,m i: m f ) = 0, (5.30) 

I. (V \ ( m i m Sx2 m ) \ 2E + P'z f -P'z? D ( 1 \ / co1 x 

b 13 (E,mi,m f ) = -( ^-Y — J —- ; — — — Re[ — : T ), (5.31) 

b 2 3(E,mi,m f ) = -b 13 (E,mi,m f ), (5.32) 

, / ^ x ,1x1(01* x9 rn 2 f 1 

b u (E,m i: m f ) = +' - 1 



Ml ' M w \E + p'J\ 2 
/^•^ l" r -> r '\ (5.33) 

^ Pz Pz ' 

b 2 A{E,mi,mf) = -bu(E,im,ixif), (5.34) 
b3li E, mi , m/) = ,!^)^_JL ¥/m ( p -- p ;.). (5.35) 

For each j, k pair, {Ij(M), I k (M)} can be seen as the coordinates of the point M in 
the (j, fc)-plane. Then the points (m u ,m c ,m t ) span a triangle in this plane whose area 
1/2 Yli,i'=i+i (mod 3) Sjk{Mi, M11) exhibits again the GIM mechanism for internal masses. Fur- 



thermore, because Im 



{K^K u ){K hV K Vt 



f) 



is antisymmetric in the exchange of i and /, 



the sum over external flavours is proportional to 

K (b jk (E,m u m f ) -b jk (E,m f ,mi)), (5.36) 

i,/=i+l (mod 3) 

which shows explicitly the GIM mechanism operating on external quarks. This pattern 
should be present in a complete computation, with pertinent modifications in the functions 
Sj k i bj k . 

All bjk but bu require the existence of the CP-even phase in Q found above, and depend 
as well on the complex reflection behaviour of the initial and final quarks, contained in p' z l 
and p'J . byi depends only on the latter, and could be 7^ only when m^rxif and E are not 
neglected in front of M w . 

Note that all bjk vanish when either or my go to zero, respecting the second point 
of the general behaviour announced at the beginning of this section. In each sector of 
quark charges (—1/3 or +2/3) the heavier masses will dominate the effect, as expected from 
intuition. Note that a two-threshold structure is present, corresponding to rrii, rrif. 

In our numerical results we use the exact values of the fonctions Sj k (M, M'). It is 
instructive, though, to show a fit for the particular case M,M' << Mw, an appropriate 
expansion for all quarks but the top, 
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and the ratio of the remaining Sj t k(M, M') to the result in eq. (\ 5.37|) is given by Si^/Si^ — > 
-J, S 1A /S lt2 -> +1.5, S 2j3 /S lj2 -> -55, S 2>4 /S 1>2 -> -J, and S 3>i / S lfl -> +|, in the same 
limit. This behaviour is consistent with the third point of the general expectations developed 
at the beginning of this section. 



5.3 Transition amplitude in terms of wave functions 



The computation of the amplitude described in the previous subsection was obtained by 
projecting over the massless poles of the propagator in the presence of the wall, once a 
flavour- changing electroweak loop was inserted. 

The result should be equivalent to the more direct one in terms of incoming and outgoing 
wave functions (wave packets), as given by eq. ( | 3.32j ), with no need to consider the full 
propagator in the presence of the wall derived in Section [|, as only the usual propagators 
were used inside loops. To prove that such is the case, consider the incoming wave function 
^inc j n e q ( | 3,24] ). The following identities! 
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lead to the equivalent expression 
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In the neighbourhood of the massless poles of the propagator of a given quark, 
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5.41) 
5.42) 



In eq. ( | 5.40| ) the first term corresponds to the residue at the massless pole of the 
"homogeneous" part of the propagator, while the term in S(p — q z ) corresponds to the residue 
at the same pole of the "non homogeneous" part of the propagator. For the particular choice 
of substraction (at q 2 = 0) in subsection |5J], the latter obviously gives no contribution. 
The conclusion found in the previous subsection, namely that only the homogeneous parts 
of the propagator participated in the amplitude under consideration, gets thus a simple 
interpretation in terms of the wave function contributions to < ip^^in^tp 1 ™ >■ And the 
equivalence of both methods to compute the amplitude is proved. 



'recall that p z 



'EP = qo in this case. 
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6 CP Asymmetry 



6.1 Reflection probability 

The task is to estimate the probability of reflecting a quark of flavour i into another flavour /, 
for a given flux of incoming particles (antiparticles), i.e. the reflected flux per unit incoming 
flux. Let us consider a finite box of length L z and unit section in the plane x — y, which 
incorporates the wall. For simplicity, we will ignore in the following the x and y directions. 
We normalize the wave functions ip mc and ip out in that box. Denote the matrix element of 
the interaction Hamiltonian, generated by the loop between the latter states, by 

< v^nticr >= A{z ^ z f) , (6.1) 

where A(i — > /) is dimensionless, see eqs. fl 5.251) , ( | 5.28 ). For any time interval T such that 
A(i - f)T/C t « oo, 

< C?|e~ iHI M" C >~ e - lEr A ^7* f) [ T/e dtS^f-^ ( 6.2) 

nf 1 1 Ul L z J -T/e 

whose modulus squared gives the reflection probability per particle (E = (Ef + EA/2 in this 
equation). Taking now the usual limit: 



[ T/e dtjz&f-** 2 
J-T/e 



— Teir5(£ { -£)). (6.3) 

T^oo 



For an incoming flux $ per unit time of particles with velocity v, the number of particles in 
the volume will be $uL 2 . The reflection probability per unit time and unit flux is then: 

p(i |A( 'r /)|2 v 2n6 t E f - 30- ( 6 - 4 ) 

f L>z 

pi 

In the infinite volume limit, J2 p f/L z — > J dpl/(2ir), and with / dp{S(Ef — Ei) = 1/v, it 
follows that 

P(z^/) = |i(^/)| 2 . (6.5) 

The number of particles reflected per unit time (or reflected flux) is $ r = \A(i — > /)| 2 $j, 
and thus related to the incoming flux through the previously computed \A\ 2 . 

Assume that the probability of the incoming state i in the unbroken phase is given by 
a distribution p(i), which depends in general on the flavour, chirality and energy of i. For 
instance, at finite temperature close to equilibrium, this distribution will be close to the 
Fermi-Dirac one, with proper modification of the dispersion relation between momentum 
and energy. We define useful CP asymmetries by: 

P {i){p{i ^ f) - p(i ^ p) 

a cp = — „ rA , 6.6 
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where the sum over i, f in the numerator implies the sum over the energies and flavours for 
the isosinglet quarks (or the isodoublet quarks , while the sum %' in the denominator acts 
over flavour and helicities and the factor 2 stems from summing over quarks and antiquarks. 
In the case of zero temperature, if we assume an equal incoming probability for all the 
flavours and helicities, and considering the asymmetry for one given energy (for instance the 
most favorable one) the formula (| 6.6| ) simplifies to: 

A cp (u) t =° = Y,&cp(u;)Iy , (6.7) 

hf 

where 

a r^=o _ - f) - p(* - 1) 

AcPWij ^ (b-8) 

is the contribution to the asymmetry of a given pair of flavours. Nfi denotes the numbers of 
flavours, 6. The colour factor has been obviated as it will finally cancel between nominator 
and denominator. 



6.2 Numerical results 

Consider first the reflection asymmetries in the "down" quark sector. The numerical con- 
tribution of the functions Sjk(Mi, My), which contain the dependence on the internal quark 
masses (I, V = u,c and t), is the following one 

£ S 1>2 (M t , M v ) = 3.82 KT 5 £ S 1>3 (M h M v ) = -1.84 1(T 5 

£ Si A (M h M v ) = 3.68 10- 5 £ S 2i3 (M h M v ) = -9.96 1(T 7 ( 6.9) 

E S 2A (M h M v ) = 5.74 ITT 7 E S 3A (M h M v ) = 6.81 10~ 7 , 

where the sums on I, V are constrained by the relation V = I + 1 (mod 3). 

The following values have been used for the particle masses: My/ = 80.22GeV, m u = 
5MeV, m d = lOMeV, m s = 200MeV, m c = 1.3GeV, m b = 4.7GeV, m t = 120GeV. Our 
numerical results are presented in Fig. [|. As expected the heavier quark masses dominate. 
This is well understood from the analytical behaviour eq. (| 5.35|) . The dominant contribution 



for charge —1/3 quarks is thus given by the s—b pair. The two spikes in the figures correspond 
to the thresholds for the external quark masses involved. Indeed, the analytical formulae in 
eq. (| 5.35| ) display such a structure, with for instance in the case sr — ► b^: 



>i 2 = 0, ^ = -{^f^, a tE = m s 



2 X2 



while for the complementary transition 6r — > si 



m s m b \ m b t . / m a m b \ m b 1 m 3 j. p _ 



6.10) 



2 2 
m s m b \ 9 /, ^, \ ( m * m b \ 3 + p _ ™ 

' 13 ~ " I IvT^ ) iCTi^' 614 + I ^7 J 2^m^7' at^-m fe . 



(6.11) 



10 r 



The reason not to add in the asymmetry the isosinglet and isodoublet quarks is the TCP' symmetry, 
eq. ( | 2.7[ ), which implies that their incoming probability p(i) are equal. The eventual effect of sphalerons to 
transform the CP asymmetry generated by the wall into a baryon number asymmetry, is anyhow different 
for doublets and singlets. We thus need to distinguish the asymmetry for doublets and singlets. 
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The total asymmetry, eq. (| 6.7|) , is clearly dominated by the SR—>bi contribution, with 
a maximum value A<£p°(c<j ~ m^) = 5 10 -26 around the m& threshold. 

Consider now charge —2/3 quarks. Strictly speaking, our computation is not a good 
approximation for this sector as in eq. fl 5.20|) we have neglected factors of mj/M^, m^/M^r, 
etc., which is not legitimate for the top quark. Nevertheless the type of GIM cancellations, 
which we justified in terms of general arguments in Sect. |5], should still hold. In particular, 
the analytical behaviour for the internal mass dependence (d, s, b in this case) , given in eq. 
( | 5.37|) , are specially accurate because all internal masses are really small in front of My/. 
The exact numerical result for these quantities is 

E S lt2 (M h My) = 2.37 10- 10 E S 1;3 (M h My) = -1.12 1(T 10 
J2S 1A (M h My) = 3.3310- 10 Y,S 2 ,z{M h M v ) = -7.37HT 12 (6.12) 
E S 2A (M h My) = 4.9 10- 11 E S 3A (M h My) = 3.45 10" 11 , 

where the sums now run on /, /' = d, s, b under the usual cyclicity constraint V = I + 1 
(mod 3). 

In this toy computation, the "up" type quarks contribution is expected to be be about 
the same as the "down" type one: while the external mass dependence gives a bigger factor, 
the internal GIM cancellation produces the opposite effect. Applying for the sake of the 
argument our analytical formulae to this sector, gives an external mass dependence at the 
dominant threshold ~ m^/M^, while the internal GIM cancellation goes as ~ m^/M^. The 
contributions of both types of quarks should not differ much. Using our formulae at face 
value results in a dominant cr — > ti contribution of order AqJ,°(u! ~ m t ) ~ 2.5 10~ 25 . In 
the present case b\ 2 is not zero, which reinforces the effect. 

As previously discussed, in a complete computation, an external suppression factor ~ 
{-ttt 1 ) 2 could disappear, leading to a large enhancement for certain flavours (up to 13 orders 

1VI W 

of magnitude for the d — s pair, still insufficient for the observed asymmetry), although this 
may be very optimistic as the pattern of unitarity triangles must persist. 

7 Conclusions 

Realistic models of baryogenesis imply the inclusion of finite temperature effects. We con- 
sider scenarios with a first order phase transition. It is interesting to disentangle whether 
temperature effects, other than the existence of the transition itself, are responsible for pos- 
sible modifications of standard CP arguments. We thus revert in a first step to the analysis 
of particle propagation at T = 0. 

The essential non perturbative element is the wall separating the two phases of sponta- 
neous symmetry breaking. The propagation of any particle of the SM spectrum should be 
exactly solved in its presence. And this we have done for a free fermion and a thin wall, 
leading to a new Feynman propagator which replaces and generalizes the usual one. It con- 
tains both massless and massive poles, and new CP-even phases, present both for on-shell 
and off-shell fermions. They stem from the tree-level reflection of fermions on the bubble 
wall. This propagator can be useful in physical process other than baryogenesis, where a 
first order phase transition is relevant. 
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When electroweak loops are considered, a CP-asymmetry is found in the one-loop reflec- 
tion properties of quarks and antiquarks. The one-loop self-energies cannot be completely 
renormalized away for on-shell fermions in the presence of a wall, as the latter acts as a 
source of momenta. The asymmetry is then found at order aw in amplitude and requires 
the interference of the mentioned CP-even phases with the CP-odd CKM ones. The de- 
sired effect is thus present without any consideration of temperature effects in the quark 
propagation. 

A complete one-loop SM calculation is beyond the scope of the present paper. We have 
instead performed a toy computation in which only the broken phase is considered inside 
the loops, in a Lorentz covariant way. The result is many orders of magnitude below what 
observation requires, and it can be expressed in terms of two unitarity triangles. They show 
the internal and external GIM cancellations and the interplay between CP-even and CP- 
odd phases. This pattern of triangles should survive in a complete computation, which we 
leave for later publication. We argue that suppression factors in external masses could be 
absent then, due to the breaking of Lorentz invariance in the loops analysis, although the 
quantitative enhancement is insufficient. 

With the expertise acquired in this paper, we face in ref.|14| the physical finite temper- 
ature case, where we discuss the common building blocks and new features of the scenario. 
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A Feynman boundary conditions for the fermion prop- 
agator 



1~T6D , ( r£T7D and ( fOBD is a unique 



We prove here that the propagator derived in eqs. 
solution of the inhomogeneous Dirac equation (|4.10|) . In a first step, the equivalence between 
eqs. ( | 4.16| ) and ( | 4.17| ) will be established. We then prove that the two 'ie" conventions are 
indeed equivalent, as stated in section (|[). Next we show, in the "space ie" convention, that 
the propagator satisfies the inhomogeneous Dirac equation. Finally, in the (Feynman) "time 
ie" convention, we demonstrate that the Feynman boundary conditions are indeed satisfied. 



A.l Equivalence between eqs. (4.16) and (4.17). 



Let us consider the difference between the right hand sides of eq. (| 4.17|) and eq. ( | 4.16|) 
The difference between the first lines of these equations is 



1 I I 1 \ 1 



ql - ql + ie \f fa f ) ql - q\ - ie \f - m 4 f -m / 



A.l) 
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where the Feynman convention is understood in all all , l/(^ — m)'s. 

The difference between the second and third lines of the same equations is 



1 1 

- T7 sa z s'y — + 



— m 



-sa z sj - 



A.2) 



— m 



qf and q l differ only in their z component, implying that $ — = ~lz(ql — q 1 )- Using 
this equality, together with — 1/^ = {£ — j*f)l/(f , it is straightforw to show that 
the sum of eqs. ( | A.l|) and ( | A.2|) vanishes, whence the equality of r.h.s. of ( | 4.16Q and of ( 



437|) . 



A. 2 Equivalence between the two "ie" conventions. 

The usual "time ie" simply amounts to q 2 — > q 2 + ie, while the "space ie" convention has 
been defined in Sec. [|. The difference between the two ll ie" conventions has been exhibited 



in eq. (| A.7|) . 

Consider the equality 



m 



+ 



$ (¥ ~ m ) $f — m § j 

Aplying it to eq. (| 4.171 ), we can separate the terms containing 1/$ and those containing 

— m. Let us dismiss momentarily the latter, which will be discussed later. 
It is possible to prove that the substitution of all \ factors in eq. ( |4.17| by 2iir^9(q[)5(ql 
(q{) 2 ) gives a vanishing result. Indeed, from 

tiWMqt ~ (q f z ) 2 ) = s\q \j°(l - sa z )6(qf)6(<lo ~ (4)% ( A.3) 

the third line in eq. ( | 4.17| ) gives then a vanishing coefficient, since (1 — sa z )(l + sa z ) = 0. 
Applying now in the second line the equality 

msjo 



sa z ) 



E + p' z 



sa 2 



A.4) 



it follows that 



(ql) 2 — (ql) 2 + ie "'qi — ql + ie 



[1 — sa z )s^n— = s(l — sa z ) — - f — = s(l — sa z )— f , ( A. 5) 



where we have used sq = q{ (since sq > by definition of s, and q z > due to the 9(ql) 
factor). Finally, using once more eq. ( | A.3| ), the second line gives 

0(q f z Wo - (gf) 2 ) J \ . , ( A.6) 

ql — q z + it 

. which obviously cancels the contribution from the first line. 

Let us now turn to the 1 / {ff — m) term. The equivalent of eq. (\ A.7|) is 
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+ m 



+ m 



+ 2vk{# ' + m)d{-q{)5{ql-{qf z ) 2 -m 2 ). ( A.7) 



The proof proceeds in a way similar to the preceding one but slightly more involved. It 
can be shown that 



+ m)6(-q{ )5{ql - (q[) 2 - m 2 ) = 8 \q \>y (l + 



p{a z + m7 ' 



A.8) 



where p{ = yq 2 — m 2 . The last bracket in ( | A.8| ) is twice a projector that plays the same 
role as 1 — sa z in the preceding derivation. 
Using 



p{ot z + 777/7°' 
s \Qo\ 



771570 
E + p' z 



sot. 



0. 



it follows that the contribution from the second line of (| 4.17) vanishes. From 



1 + 



p{.a z + 7717 \ 1 + set 2 



s 5b 



1 + 



ms7o 



1 + 



p{a z + 7777° 



A.9) 



A.10) 



E + p' z 

it is possible to check that the third line cancels the contribution from the first one, which 
ends the proof for the for 1 / — m) term. 

The — 777) and may be treated in a similar way leading to the anounced 

conclusion: the quark propagator in eq. ( |4.17| ) has the same value whichever "ze" convention 
is used. From the equality derived in the preceding subsection, this applies also to eqs. ( 
47T6|) and (POT). 



A. 3 The inhomogeneous Dirac equation. 

We prove here that the propagator verifies the inhomogeneous Dirac equation. Let us con- 
sider the equation ( | 4.10| ). Multiplying both sides by 770, the Dirac operator becomes 

O(-Qi0 + e(Q(ip-m). ( A.ll) 

We use the "space ie" convention, so that \ has a pole in the q{ complex plane below the 
real axis. Its Fourier transform has then a pole above the real axis, 

Jdqf z e^jj^e(-Cm-\ ( A.12) 

Similarly, l/(^ — 777) has a pole above the real axis, leading to 

/ d ^ eiqki J^ K - m l _1 - ( A.13) 

As a consequence, the operator ( | A.ll| ) applied to Tn/tf*($' — m) gives zero. It follows 
that the second and third lines in eq. ( | 4.17|) vanish under the action of the Dirac operator. 
Its action on the first line gives 
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- f ~. , . + f \ . = 2i*6(qf - ql). ( A.14) 

qi - q\ + ie qi — q\ — ie 

This ends the proof. The same argument holds when the Dirac equation is applied on 
the right hand side of the propagator. In this case, it is more convenient to perform the 
analysis on eq. ( | 4.16Q , sticking again to the "space ie" convention. 

A. 4 Feynman boundary conditions. 

In the "time ie" convention, it is trivial to proof that the propagator verifies the Feynman 
boundary conditions. Indeed, the latter state that the positive frequencies are "retarded" 
while the negative ones are "advanced" . By definition, the "time ie" prepscription achieves 
it, as it amounts to the replacement q$ — > q$ + ie, which leads to poles in the go complex 
plane below (above) the real axis when its real part is positive (negative). 
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Figure 4: Asymmetries for incoming down quarks, as a function of energy. 
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